Electron correlation resonances in the transport through a single quantum level 
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Correlation effects in the transport properties of a single 
quantum level coupled to electron reservoirs are discussed the- 
oretically using a non-equilibrium Green functions approach. 
Our method is based on the introduction of a second-order 
self-energy associated with the Coulomb interaction that con- 
sistently eliminates the pathologies found in previous pertur- 
bative calculations. We present results for the current-voltage 
characteristic illustrating the different correlation effects that 
may be found in this system, including the Kondo anomaly 
and Coulomb blockade. We finally discuss the experimental 
conditions for the simultaneous observation of these effects in 
an ultrasmall quantum dot. 

The transport properties of a quantum dot coupled to 
electron reservoirs have received considerable attention in 
recent years both from the experimental jlj and the the- 
oretical [|| sides. With the advent of nanotechnologies, 
an ideal experimental device where tunneling takes place 
through a single electronic level within the dot is now fea- 
sible. Different theoretical papers [Q-Q have suggested 
that the tunneling through this single quantum level can 
be appropriately described by the Anderson hamiltonian 

H = Y^ ^n Qa + T kKl c o<? + C L C D 

+ Y,^ + ^Xa + Un a] n all (1) 

where eo represents the quantum dot level, U is the effec- 
tive Coulomb repulsion within the dot, e£ + fx u with v = 
Left or Right denotes the reservoir single-particle ener- 
gies, fi L — n R = eV being the applied bias and T£ the 
coupling between the reservoir states and the quantum 
dot level. 

Hamiltonian (1) is expected to describe the main cor- 
relation effects associated with the quantum dot. In par- 
ticular, a Kondo- like resonance near the Fermi energy 
should be reflected in the current-voltage characteristic 
at low voltages, when the dot level is nearly half-filled. 
On the other hand, a Coulomb blockade effect should 
lead to a reduction in the conductance up to voltages 
large enough to overcome the Coulomb repulsion in the 
dot level. Previous theoretical efforts have mainly been 
concerned with one of these two effects. In references 
the Kondo-effect is studied using a variety of tech- 
niques both in the linear and non-linear regimes, whereas 
the charging effects has been analyzed in references . 
The aim of this paper is to present a solution for this 
model which would allow us to describe accurately all 



the physically relevant regimes, showing how many-body 
resonances appear in the current-voltage characteristic. 
We pay particular attention to the situation where fea- 
tures related to the Kondo and charging effects may be 
simultaneously found. Based on these results, we shall 
also discuss in detail the experimental conditions neces- 
sary to observe those resonances in semiconductor nanos- 
tructures. 

The transport properties of this model can be analyzed 
by using the Keldysh formalism where the retarded, 
G r , and the distribution, G H , Green functions are de- 
fined as follows 

Gii( w ) =-* [ < 4Wc*(0) + Ci(0)cUt) > e lwt dt 



(2) 



G tji w ) = i I < 4(t)ci(0) > e lwt dt. 



These Green functions can be calculated by starting 
with the L and R reservoirs decoupled from the dot 
level. This case defines the unperturbed Green functions 
G {0 \w). ThenG(w) can be obtained by coupling the dot 
level to the reservoirs and by introducing the self-energies 
H r and X H which take into account the Coulomb cor- 
relations within the dot. 

Once the different Green functions are obtained from 
the corresponding Dyson equations ||, the current in- 
tensity I v between reservoir v and the dot level, is given 
by 
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The crucial point in order to solve this problem is to 
find a reasonable approximation for the self-energies. For 
U = 0, S = and the exact G(w) can be easily obtained 
using conventional Green function techniques [[10[ . 

The effect of a finite U can be included by using per- 
turbation theory in U. This perturbative approach has 
been extensively analyzed by Yamada and Yoshida |ll]] 
and Zlatic and Horvatic |L2j (hereafter referred to as YY 
and ZH) for the equilibrium Anderson Model (zero ap- 
plied bias). Hershfield et al. [Q have extended this ap- 
proach to the noncquilibrium case by calculating Y7 Qty (w) 
and T,q~(w) up to second order in U. 

(2) 

In YY and ZH approach, T, J(w) is calculated from 
the second-order diagram shown in figure 1 (inset), where 
each Green function line is a Hartree Fock dressed prop- 
agator, whose retarded part is given by: 
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(4) 



where e 0a = e + U < n QS > and T u = 2ir I T k I 2 
5(w - e k ). 

Although this approximation gives a good descrip- 
tion of the electron correlation effects in the symmet- 
ric case (eo = —17/2), it presents some drawbacks when 
one moves away from this condition. In the equilib- 
rium case this is clearly illustrated by the failure of 
the Friedel-Langreth jl3| (hereafter referred as FL) sum 
rule, that relates the "impurity" charge < riQ a > to the 
phase shift created by the impurity at the Fermi energy 
r)(0) = Im[lnG r 0a (0)}. 

Figure la shows the degree of fulfillment of the FL 
sum rule, < no a >= — hvifyi using the second order self- 
energy £g^(u>) in the YY-ZH approach. In this figure 
we compare < no a > with ——r](0) as a function of the 
dot level eo, for U/T — 2Air, neglecting the frequency 
dependence in r„, and taking = Tr = T/2. 

(2) 

Improvements over the above T, Q J are not easy to ob- 
tain. For instance, one could insert the full self-consistent 
dressed propagators instead of the Hartree-Fock ones for 
calculating the diagram in figure 1. This sort of self- 
consistent perturbation theory is charge conserving fli] ] 
and would verify the FL sum rule. However, it has been 
shown that it leads to a poorer description of the quasi- 
particle spectral-density [FBI . 

Our proposal to improve the YY and ZH approach is 
the following: instead of using the Hartree-Fock solution 
as the initial one-electron problem for the calculation of 

(2) 

S 0o . , we use a different self-consistent field as a starting 
point. In this initial one-electron problem the Hartree- 
Fock level, eoo-, in Eq.(4) is replaced by an effective level 
£e//cr, which is determined by imposing self-consistency 
between the initial and the final dot level occupancies. 
As shown in figure lb, this way of calculating the second 
order self-energy gives a solution that closely verifies the 
FL sum rule. 

Notice that this procedure is formally correct within 
perturbation theory: we have decomposed H — H e ff + 
V, where H e ff is the one-electron part of hamiltonian 
(1) with eo replaced by e e // CT , and we treat V = H — 
H e jf as a perturbation. The fulfillment of the FL sum 
rule and the condition of charge consistency provides an 
appealing interpretation for e e ffa- it must behave like 
eo+ ReT,Q a (0), the effective potential at the Fermi energy. 
Obviously, for the symmetric case our approach coincides 
with that of YY and ZH. 

This procedure allows us to get a good description 
of the Kondo-like peak appearing around the Fermi en- 
ergy and the overall behavior of the spectral density for 
a broad range of parameters (eo, U/T). However, for 
high values of U/T and far from the symmetric situation, 
the position and the spectral weights of the resonances 
around eo and eo 4- U are not so well described 119] . The 



reason is that does not yield the atomic limit when 

w ~ U 3> r. To improve the solution given above in this 
region we introduce the following self-energy |ll| 
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with 



a = 



(1- < n g >)U + ep - e e // 
< n og . > (1- < n m >)U 2 ' 



which has the virtue that it yields the appropriate atomic 
limit ]l6| , |r7t , for w ~ U and T/U — > and behaves like 



,r(2) 



for U/T -> 0. 



Let us now turn our attention to the non-equilibrium 
situation. For the sake of simplicity, we shall assume 
that in Eq. (1) [i L — —[i R — eV/2. Following our ap- 
proach for the equilibrium Anderson model, we obtain 
the different nonequilibrium self-energies, Sq ct and S^T, 

R 
l eff 



by introducing the effective levels e e // CT , Me// an d fi^ 



in the initial one-electron hamiltonian 



H 



(transport) 
eff 



E 



u.k.a 



(6) 



As a natural extension of our procedure for the equi- 
librium case, e e ff<T) Heff an d A^ff/ are determined by 
imposing self-consistency in the dot level charge and in 
the currents II and Ir defined by Eq. (3). 

Eq. (6) allows us to calculate the one electron Green 
functions G e/ 1 that will be used to obtain zZ QtT and 
S 0(T by means of a second order perturbative calcula- 
tion, and finally Y,Q a (w) and Tiq~(w) by using the equiv- 
alent of Eq. (5) for the non-equilibrium case. 

Let us remark that our method eliminates, in a natural 
way, the pathologies that arise when the perturbation is 
performed upon the Hartree-Fock solution. In the non- 
equilibrium case these pathologies are reflected in the 
unphysical behavior of the dot level charge as a function 
of eo (around eo ~ —U/2); the currents 7l and Jr exhibit 
also unphysical features Q] including violation of current 
conservation (II ^ Ir). A detailed comparison of the 
results given by the different approaches, together with 
a more comprehensive discussion of our method will be 
presented in a future publication. 

In order to illustrate the kind of results obtained by 
our method we have considered two different situations, 
with the dot level below (Fig. 2) and above (Fig. 3) 
the Fermi energy respectively. In the second case the dot 
level is almost empty at zero bias. 

Figure 2 shows the current intensity, J, and the dif- 
ferential conductance g = dl/dV, as a function of the 
applied bias, for eo = — 7r and U = 0, 7r/2, tt, 2tt, and 



47T, where all the energies hereafter are measured in units 
of r. Our results for g show a single broad peak around 

V = 2tt for U = 0; this corresponds to the dot level 
€q = — 7T crossing the right reservoir chemical potential 
fi R = —V/2. For small U (U < |eo|) this one-electron like 
resonance is shifted and adopts a somewhat asymmetric 
shape. On the other hand, for large U (see case U = 4-7T in 
Fig. 2) three different features are clearly present in the 
current-voltage characteristic. The conductance peak at 

V = is related to the Kondo resonance appearing in the 
dot spectral density around the Fermi energy j|, while 
the peaks at V = 2tt and V = 4ir correspond to the 
crossing of the dot "ionization" and "affinity" levels, eo 
and eo + U, with the reservoirs chemical potentials. We 
should comment that these two peaks reflect the charg- 
ing effects associated with the dot level: the second level 
can only be filled when the applied bias overcomes the 
repulsion between the second and the first electron. For 
the particular case U = 2tt we have a symmetric prob- 
lem with the two dot levels crossing the left and right 
chemical potentials at the same bias, V — 2ir, leading to 
a single "charging effect" peak. Notice that in this case 
the conductance at V = reaches its maximum value 
2e 2 /h. Finally, the case U = it illustrates the transi- 
tion between small and large U , with a new correlation 
structure arising due to the overlap between both Kondo 
and charging effects, the resonances merging into a broad 
peak around V = 0. This transition is also apparent in 
the current intensity I(V) shown in Fig. 2, where the case 
U = 7T defines the border between the highly correlated 
limit (Kondo and charging effects completely separated) 
and the one-electron like behavior. 

Figure 3 shows I(V) and g as a function of V for 
eo = 37T and U = 0, 2ir, 4ir and 6ir. In these cases 
the dot level is above the Fermi energy, and for V = 
no Kondo-like peak appears in the spectral density of 
states. As before, the two peaks appearing in the differ- 
ential conductance for U > 2tt are related to the filling of 
the ionization, eo, and affinity, eo + U, levels as a function 
of the applied voltage. As shown in figure 3b, the voltage 
difference between the two peaks for U > 2tt is approxi- 
mately equal to 2/7. On the other hand, the conductance 
peak at V = 2e is reduced to nearly half its U = value 
when the interaction increases. It is also worth noticing 
that the I — V curve in this case exhibit a step-like be- 
havior which is more pronounced for increasing values of 
U. The relative height of these steps is in agreement with 
calculations based on a simple atomic-like model for the 
dot Green functions ||. 

Figure 2 and 3 show the different kind of current- 
voltage characteristics one gets for the single dot level 
described by hamiltonian (1). For U sufficiently large, 
two resonances associated with the ionization and affinity 
levels can always be observed in the differential conduc- 
tance. In this case correlation effects are very important 
giving rise to an additional resonance at V — when the 



ionization level is initially filled. 

The results found in this paper suggest that the Kondo- 
like resonance and the peaks associated with the charging 
effects of the quantum level could be found simultane- 
ously in a quantum dot, provided that the temperature 
is sufficiently below the Kondo temperature. As an ex- 
ample, we consider a GaAs-dot of size L ~ lOOA both 
in the vertical and lateral directions, sandwiched between 
two GaAs- wires and two Al Ga As-barriers, as those stud- 
ied in reference . A single bound level around 50meV 
is found by solving numerically a simple double barrier 
model with a barrier height ~ 300 meV . The parameter 
U can be evaluated as the Coulomb integral for the wave 
function corresponding to this bound level, which roughly 
yields U ~ e 2 /eL ~ 15meV, in agreement with the dot 
classical capacitance energy [|| . Then, varying the AsGa- 
wires doping around 10 18 /cm 3 , one could get an experi- 
mental device close to some of the theoretical cases ana- 
lyzed in figures 2 and 3. In particular, for n = 10 18 /cm 3 
and a barrier width of 30 A, we find Cq/T ~ —3 and 
U/T ~ 8, not far from one of the cases presented in 
figure 2. We should recall that finite conduction band 
effects are not included in our present calculation; one 
should always keep in mind that for V sufficiently large 
the effect of the bottom edge of the semiconducting wires 
would be present in the I — V characteristic. 

In conclusion, we have presented an accurate solution 
for the many-body problem of a single dot level between 
two biased reservoirs. Our results show the experimen- 
tal conditions one should achieve to observe in the dif- 
ferential conductance a peak related to the Kondo-like 
structure in the density of states, and a complementary 
structure associated with the charging effects of the dot 
level. By adjusting appropriately the dimensions and the 
doping of a quantum semiconducting dot, we have shown 
how the different peaks could be determined by measur- 
ing the differential conductance. 
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FIG. 1. Fulfillment of the Friedel-Langreth sum rule as 
a function of the dot level position using the second order 
self-energy as calculated: (a) starting form the Hartree-Fock 
solution and (b) imposing self-consistency in the dot level 
charge. The full line corresponds to the dot level charge 
< nocr > and the dotted line corresponds to — ^??(0). Inset: 
second order diagram used to calculate the self-energies. 

FIG. 2. Current I and differential conductance g obtained 
using our method for eo = — w and (7 = (a), -k/2 (b), ir (c), 
2ir (d), and 4ir (e). All the energies are measured in units of 
the elastic decay rate at the dot T. 



FIG. 3. Same as figure 2 for eo = 37r and U = (a), 2tt 
(b), 4tt (c), and 6tt (d). 



